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^ Abstract 

In this paper, the Brownian motion theory has been applied to the turbulent relative 
particle dispersion problem. The fluctuating pressure forces acting on a fluid particle are 
taken to follow an Uhlenbeck-Ornstein process and their correlation time appears to be 
very plausibly identifiable to go inversely proportional to the flow Reynolds number R e . 
This connection leads to the result that the Richardson- Obukhov scaling holds only in the 
infinite R e limit and disappears otherwise. Further, the Lin-Reid [TT] conjecture regarding 
the connection between the pressure force fluctuation parameter and the energy dissipation 
rate in turbulence is confirmed and a -R e -dependent explicit relation between the two is 
determined. 
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1. Introduction 



Richardson [T] proposed that turbulent diffusion should be characterized by the distance 
between neighboring particles, because if the interparticle distance is within the inertial 
range, one may expect to find universal super diffusive behavior in the relative particle dis- 
persion process, which may be interpreted in terms of an interparticle separation depen- 
dent turbulent diffusivity. From a purely empirical analysis of atmospheric dispersion data, 
Richardson pQ then showed that the turbulent diffusivity defined by the rate of increase of 
the mean square interparticle separation distance, 

wo a) 

goes like 4/3 power of this distance, 

D ~ ([R(t)} 2 ) 2/3 . (2) 

Obukhov [2] showed that Richardson's relation (2) can be derived via Kolmogorov's [3] 
theory for homogeneous isotropic 3D fully developed turbulence (FDT). When the interpar- 
ticle separation is within the inertial range, Obukhov |2j gave 

D~ei(iR(t)] 2 f 3 (3) 

e being the mean energy dissipation rate, hence the name Richardson- Obukhov (RO) scaling. 
On the other hand, Richardson's pQ formulation connects with the universal aspects of FDT 
actually stronger than Kolmogorov's formulation since the unphysical effects due to sweeping 
by larger scales are precluded from the onset in Richardson's pQ formulation. 

The RO scaling result has received, however, little experimental support due to the diffi- 
culty of performing Lagrangian measurements over a broad enough range of time and with 
sufficient accuracy. Even in recent laboratory experiments (Ott and Mann |3], Sawford [5], 
Bourgoin et al. |6j, Salazar and Collins [7]) with high-speed photography to track particles 
and in numerical simulations (Yeung [8], Boffetta and Sokolov [9]) with the highest possible 
resolution for homogeneous isotropic turbulence, it is known to be hard to obtain an ex- 
tended range with the RO scaling. Indeed, Bourgoin et al. [5] reported that the RO scaling 
may not be observable in laboratory experiments even at high Reynolds numbers unless the 
initial particle separation is significantly small (so as to preclude the ballistic regime that 
currently seems to dominate the observable regime). The difficulty of achieving the RO 
scaling in laboratory experiments and numerical simulations appears to be due to, 

* contamination of the inertial range by dissipative effects at the ultraviolet end and by 
the external forcing effects at the infrared end of the spectrum caused by inadequate 
scale separation; 

* persistent memory of initial separation; 

and is therefore basically traceable to the finiteness of the Reynolds number as well as the 
observational domain. The purpose of this paper is hence to shed some light on this aspect by 
using the Brownian motion theory applied to turbulent relative particle dispersion problem 
(Obukhov [10]). 
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2. Application of the Brownian Motion Theory 



A fluid particle in a turbulent flow moves in response to fluctuations in the pressure of 
the surrounding fluid. The assumption of short-range correlation of the fluctuating pressure 
forces in a turbulent system leads to the applicability of the Brownian motion theory to 
the turbulent dispersion problem (Obukhov [10]). Obukhov [10] proposed that the motion 
of a fluid particle mimics a Markov process and is therefore described by a Fokker-Planck 
equation. The fluid particle is assumed to be subjected to successive small impulsive forces, 
like those experienced by a Brownian particle. These assumptions also hold for the relative 
motion of two particles in view of the negligible effect of large-scale motions in this situation 
(Lin and Reid [II]). 

Consider the relative motion of two particles released at the same point with the same 
velocity in an unbounded fluid in a state of stationary homogeneous turbulence. The particles 
then drift apart due to different initial accelerations. We then consider the statistical average 
over an ensemble of such pairs. This problem is governed by the following initial-value 
problem, 



^ + fiV = a (() (2) 



with, 



t = : R = 0,V = 0. (3) 

Here, /3 is a constant coefficient of resistance and a (t) is the fluctuating pressure force which 
follows a stationary, random Markov process 
The solution of (l)-(3) for a particle pair is 

V{t)= /V^a(£)de. (4) 
Jo 

The mean square of relative velocity of an ensemble of such particle pairs is 

< [V (t)f > = e~ m f drj f di e^ +ri) < a (f) a (77) > . (5) 



jo jo 

We now assume that a (t) has an auto-correlation function given by (Uhlenbeck and 
Ornstein p.J), 

<a(t') a (t")>=a 2 e-^ t '- t "^ (6) 
A being the inverse of the correlation time. Large A corresponds to white nois <E(the 

case 

considered by Lin and Reid [H]) while small A corresponds to persistent pressure fluctuation 
correlation. 



2 The white noise assumption for the fluctuating pressure forces in a real turbulent flow situation does not really seem to be 
justifiable. 
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Using (6), (5) becomes 



<[V(t)] 2 >- 



(j- 

T 



A 



-2/8* 



f3\\ + f3 X-f3J\ 2 -f3 



+ 



2A 



,-G8+A)t 



(7) 



For large A, (7) becomes (Lin and Reid [TT]). 

a 2 /A 



<[v(t)Y> 



(3 



(1 - e" 2 ^) 



which leads to 



2 I y ) / . )'/ < I 



(oVA) 



(9a, b) 



/3t =>- oo. 



(8) or (9) shows that, in the absence of dissipation ((3 = 0), the particle pair experiences 
a runaway drifting motion due to the cumulative effect of the fluctuating pressure forces 
acting on the particle pair which follow a stationary process. 

For small A, (7) becomes 



<[V(t)] 2 > - (j ) j 2 (1 + - 2e~») 



(10) 



which leads to, 



a 2 t 2 ,/3t < 1 



<[V(t)] 2 > 



a- 

J 2 



(11a, b) 



(3t =>• oo. 



Comparison of (11) with (9) (see Figure 1) shows that the mean square of relative velocity 
initially grows more rapidly for the short-range correlation case than for the long-range 
correlation case. 

Next, using (4), (1) and (3) give 



R(t) 



-/3(*-0 



The relative particle dispersion is then given by 



(12) 



< [R (t)] 2 > = ^e~ 2 ^ f d V f di < a (0 a fa) > 

P Jo Jo 



+ wl dr] I d t <a ® a W > ~i e ' m J dr] J d ^ eK <a (o«(^) > • 

(13) 

Using (6), (13) becomes 
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<[V(t)] 2 > 




Figure 1 . Effect of the fluctuating pressure force correlation 
time on the mean square velicity of the Brownian 
particle. 
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< [R(t)} 2 > 



(oVA) 

/3 2 



A / 1 



-2/3t 



P\\ + P X-f3\ 2 -f3 



+ 



2A 



-(A+/3)t 



+ 2 



>7A) 

/3 2 



1 1 



-At 



(oVA) 

/3 2 



2A + /3 
L/3(A + /3) ' (A + /3) 



+ 



-(A+/3)t 



(14) 



(A-/3) /3(A-/3) 



For large A, (14) becomes (Lin and Reid [IT]). 



<[fl(*)r> 



2. . ^ 2 A 



/3 3 



(l- e - 2 ^)+2^-4^ (1-e-^) 







0» 



(15) 



which leads to 



<[i?(t)] 2 > 



3 V A 



;i6a, b) 



t, /3£ =>■ oo. 



I /3 2 

(16a) represents the Richardson [1] scaling and signifies the validity of a Brownian motion 
theoretical framework for the relative particle dispersion problem. Indeed, if one defines a 
relative diffusivity D by 



(17) 



then (16a) gives (Lin and Reid jTTJ) 

2/3 / 2 \ 1/3 



c= (x) t2= G) (x) (< wi 2> ) 



2/3 



which is Richardson's [TJ other scaling result implying that the relative diffusivity D increases 
as the 4/3 power of the particle separation distance. 

In view of the consistency of Richardson [JJ scaling with Kolmogorov [3] scaling (which 
is valid in the infinite flow Reynold's number R e limit), it therefore appears very plausible 
to take the Brownian parameter A to be proportional to the flow Reynolds number R e , 

A = kR e (19) 

k being a constant. It is of much interest to note that the connection between the white noise 
assumption for the fluctuating pressure forces and the Kolmogorov [3] large R e assumption 
indicated by (19) was querried earlier by Corrsin 



G 



For small A, (14) becomes 



2e~ pt ) 



(20) 



which leads to 



<{R{t)] 2 > 



(^/A) V;J 



4 

> 2 A) 



At 4 , pt < 1 
t, fit =>- oo. 



(21a, b) 



(21a) shows that the Richardson pQ scaling holds only in the infinite R e limit and disappears 
at finite i? e 's. The relative diffusivity D, defined as per (17), is now given by 



D 



1 (a 



2 \ A 



At 3 



4A 



2\ 1 1/4 



([R(t)] 2 >y 



(22) 



which shows that the relativity diffusivity D, in the small R e limit, increases as the 3/2 
power of the particle separation distance. 



3. Relation between the Brownian Motion Parameters and the 
Kolmogorov Parameters 

In the Brownian motion framework, as indiated by (14), the relative particle dispersion 
depends on the parameter (a 2 / A). As pointed out by Lin and Reid [TT], the parameter (cr 2 / A) 
looks very much like the energy dissipation rate parameter e in the Eulerian Kolmogorov 
[3] framework {a 2 /\ and e even have the same dimensions). Lin and Reid [11] therefore 
went on to contemplate the possibility that " (a 2 / A and e) might be identifiable, apart from 

a 2 / A 

a constant factor. But, in general, one can only surmize that perhaps = / (R e ) " . We 

will now show that (19) indeed underscores such a possibility and hence contributes toward 
providing a Brownian motion basis for the Kolmogorov [3] theory. 

Multiplying equation (2) by V and taking the average over an ensemble of such particle 
pairs, we obtain 

\j t <[V{t)f> + P<[V{t)f> = <a(t)V(t)>. (23) 
Upon substituting (4) and using (6), we obtain 

\± <{V{t) f > + S< { V {t) f> (24) 
Further, on assuming a stationary state in the limit t — > oo, (24) yields 

e ~^mi 2> =ttS) (25) 
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On using (19), (25) leads to 



(^)=/ ( /y~i + *(#) (26) 



^Re 

which confirms the Lin-Reid [11] conjecture and provides one with an explicit determination 
of / (R e ) hinted at by Lin and Reid [TT] . 



4. Discussion 



There is speculation that the difficulty in obtaining an extended range with RO scaling 
in both laboratory experiments and numerical simulations is due to the finiteness of the flow 
Reynolds number R e in these situations. In this paper, the Brownian motion theory hs been 
applied to the turbulent relative particle dispersion problem to shed some light on this issue. 
The fluctuating pressure forces acting on a fluid are taken to follow an Uhlenbeck-Ornstein 
process and their correlation time appears to be very plausibly inversely proportional to 
R e . This confirms the connection querried earlier by Corrsin [T3] between the white noise 
assumption for the fluctuating pressure forces and the large R e assumption in the Kolmogorov 
[3] theory. This connection leads to a totally reasonable result that the RO scaling holds 
only in the infinite R e limit and disappears otherwise. Further, the Lin-Reid [11] conjecture 
regarding the connection between the pressure force fluctuation parameter and the energy 
dissipation rate in FDT is confirmed and a i? e -dependent explicit relation between the two 
is determined. 
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